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Abstract. Complex period variations of five W UMa type binaries (AB And, 00 Aql, DK Cyg, V566 Oph, U Peg) were 
^ ' investigated by analyzing their 0-C diagrams, and several common features were found. Four of the five systems show secular 

period variations at a constant rate on the order of IP^cc/Pl ~ 10 ' y '.In the case of AB And, 00 Aql, and U Peg a high- 
amplitude, nearly one-century long quasi-sinusoidal pattern was also found. It might be explained as light-time effect, or by 
some magnetic phenomena, although the mathematical, and consequently the physical, parameters of these fits are very prob- 
lematic, as the obtained periods are very close to the length of the total data range. The most interesting feature of the studied 
O-C diagrams is a low amplitude (~ 2 - 4 x 10"^ d) modulation with a period around 18-20 yr in four of the five cases. This 
\ phenomenon might be indirect evidence of some magnetic cycle in late-type overcontact binaries as an analog to the observed 

activity cycles in RS CVn systems. 
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1 . Introduction period variations. The latter may also be of various types, e. g. 

caused by the evolution of the system, mechanical and/or ther- 

The moments of minima are one of the most fundamental ob- , • i a- , a ^ • ui . 

modynamical enects, or due to variable magnetic activity. Of 

, servables of eclipsing binary systems and, due to their rela- ^^^^^^^ ^ ^^^^^^ ^-^^^^^^^ ^^^^^^^ ^^^^^ 
^ . lively easy measurement, a huge pile of data of this type was ^^^^ ^.^^^ generating a very complex variation in the orbital 
" " ' collected by observers during the last century. However, so far pgj-jQjj 
there is no a commonly accepted and straightforward interpre- 
tation of the O-C diagrams constructed from them. Additional I" respect to the involved time scales we can distinguish 
difficulties arise when one tries to disentangle and identify the *e long-term variations (on a nuclear or thermal time scale) 
physical mechanisms causing the observed period variations al- *e short-term variations characterised by a decade-long 
though this would be the aim. For example, in spite of the fact cycle length. In our recent study we focus on the short-term 
that the period variations are common features of contact bi- orbital period variations of several W UMa systems, because 
nary stars (Kreiner 1977; Maceroni & van't Veer 1996), their *e real O-C observations can reveal only these effects and the 
origin is not well understood yet. long-term ones require a different (statistical) approach. 

These observed variations can be classified into two sub- 
sets: (i) in the first part are the apparent period variations due 
to light-time effect (LITE) caused by a distant third or further 

body', and (ii) the second part comprises the inherent physical r, »i • » i » » j w t » » 

-' ' ^ ^ ^ ^ ^ -' Recently, the internal structure and evolution of contact 

Send offprint requests to: T. Borkovits, W UMa systems are not well-known or else are poorly under- 

e-mail- borko@alcyone bajaobs hu stood. The evolution theories of these binary stars yield secular 

' The other fundamental class of the apparent period variations, as well as shorter time scale quasi-cyclical orbital period vari- 

namely, the effect of the apsidal motion (AM) is not considered here ations. Although these theories suffer from a lot of problems 

due to the circular orbit of the systems. their predictions can be used as guidelines (Kreiner et al. 2003; 



Long-term effects 
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Kahler 2004) to better utiderstatid the behaviour of these sys- 
tems. 

All the currently available theories predict a secular or- 
bital period increase of such systems caused by mass trans- 
fer between the components. Some theories, e. g. the Thermal 
Relaxation Oscillation (TRO) (Lucy 1976; Webbink 1976) or 
discontinuity (DSC) theory (Lubow & Shu 1977), give rise 
to different periodic and secular variations of the orbital pe- 
riod and of the observable light curves. However, Mochnacki 
(1981); Kahler (1986), as well as Hazlehurst (2001) showed 
that none of these theories agree with the observations or seem 
to violate the second law of thermodynamics. The Angular 
Momentum Loss (AML) theories are based on the assump- 
tion of the angular momentum loss via magnetic breaking on 
stellar winds (van't Veer 1979) or on the magnetic fields of 
the companion star (Mochnacki 1981), or even on the simple 
mass and angular momentum loss via stellar wind (van't Veer 
&Maceroni 1989). 

Kahler (2004) integrated the evolutionary equation of a 
model contact system and found that the systems are oscillat- 
ing both with and without losing their contact configuration. 
During these oscillations the period shows (quasi-)periodic 
variations on a time scale of several million years. The origin 
of this period variation is once again mass transfer between 
the components. Note that Qian (2001a) suggested that sys- 
tems with q < 0.4 show secular period decrease and systems 
with q > 0.4 secular period increase (q is the mass ratio), an 
efi'ect that - if real - has not yet been interpreted theoretically. 

Independent of any uncertainties in them, all of these the- 
ories predict the period variations on a time-scale comparable 
to the nuclear evolution time- scale of the primary components 
i. e. some billion years which is far from what is observed in 
the O-C diagrams. 

Short-term effects 

Short-term effects are the most exciting ones for us because 
they can be studied through the observed times of minima and 
the O-C diagrams. 

The LITE that cause apparent orbital period variations was 
studied in many contact binary systems, and one of the most 
recent extended studies was pubUshed in Borkovits & Hegediis 
(1996), also includes 18 Algol and contact systems. They found 
that only four of the 18 cases studied might have a third (or 
fourth) companion with high probability. Remarkably two of 
the studied five W UMa systems, namely AB And and AK Her, 
could have a third component compared to the two Algol type 
systems from the fifteen ones considered. 

However, Hendry & Mochnacki (1998) found that AK Her 
does not show spectroscopic evidence for a third companion in 
the system. Note that this star has a wider, visual companion 
at a separation of p = 4'.'7, nevertheless this companion cannot 
be responsible for the decades-scale cycUc variation of the O- 
C. Since the time coverage of their observations was only a 
few days, the gamma velocity was not changed, while the lack 
of the third body's spectroscopic line can be explained by a 
low luminosity object (e. g. a white dwarf). They also found 



from spectroscopy that SW Lac and V502 Oph have a third 
body, although the O-C diagrams did not show their presence 
clearly. These facts might be explained by a high inclination 
effect. This contradiction between the O-C diagrams and the 
spectroscopic observations has not been solved yet. 

Note that the light-time effect explanations of the O-C di- 
agrams are uncertain in those cases when less than one revo- 
lution of the third body could be observed till now. One im- 
portant example for this latter statement can be the case of AK 
Herculis (Rovithis-Livaniou et al. 1999). At the same time, in 
the case of VW Cephei the period of the LITE orbit deduced 
from the O-C diagram is in very good agreement with the as- 
trometric observations (Hershey 1975; Heintz 1993), although 
the amplitude of the LITE and, consequently, the projected 
semi-major axis of the third body orbit was found to be incon- 
sistent with these other measurements (see Kaszas et al. 1998). 
Nevertheless, independent of the possible reasons for this in- 
consistency, it seems evident that at least a significant part of 
the cyclic O-C variation arises from the astrometrically, and in- 
terferometricaUy identified tertiary. Although it is clear that the 
LITE studies should be improved, it is not an easy task consid- 
ering that the typical time-scale of the revolution period of the 
third body is as long as several decades or centuries. Thus, it 
should be emphasised that independent spectroscopic affirma- 
tions are inevitable before one can approve any LITE solution 
constructed from O-C diagrams. 

A further reason for period variations on a short time-scale 
can be the magnetic activity of the contact binary system. There 
are significant pieces of evidence for their magnetic activity, 
such as 

(i) Doppler-tomography (Barnes et al. 2004; Hendry & 
Mochnacki 2000; Maceroni et al. 1994); 

(ii) O'Coimell-effect and their spot-explanation in many sys- 
tems; 

(iii) their X-ray emission deviations (Stepien et al. 2001); 

(iv) Ha emission in some W UMa systems (e. g. Hrivnak et al. 
1995; Csizmadia et al. 2004). 

From an extended study of the period variations of a sample 
containing roughly 100 Algol systems. Hall (1989) concluded 
that the spectral type of systems showing cyclic period vari- 
ations is later than F5. This earlier hypothesis was supported 
during the last decade by several investigations. 

Magnetic activity is a manifestation of strong magnetic 
fields present in the system, which interact with the matter 
of stars in plasma state and could change the inertia tensor of 
component(s) through this interaction. Since the inertia tensor 
is a part of the total angular momentum required to be con- 
stant, the variation in the magnetic field structure might in- 
fluence the orbital period of the system in this way. This ef- 
fect is known as "the interchange of the magnetic and kinetic 
energy" (Applegate 1992; Lanza & Rodono 2002). The time- 
scale of this so-called Applegate-mechanism is several years or 
decades, and its amplitude may be as large as 0.006 days. This 
value is about 10-20 times higher than the accuracy of one min- 
imum time observation. 

Other reasons for period variations can also be imagined, 
e. g. sudden mass exchange via large flares, gaining a planet. 
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etc. Time-dependent variable long-scale current in the common 
convective envelope could also cause variable orbital periods, 
while all of these effects can affect the O-C diagram. 

Since now more minima observations are available than a 
decade ago and further theoretical calculations on the perturba- 
tional effects of a third body were carried out (e. g. Borkovits 
et al. 2003), we decided to improve and re-analyse the O-C 
diagrams of some eclipsing contact binary stars. In Sect. 2 we 
describe the method and in Sect. 3 we study these systems in- 
dividually. The variations found in the O-C diagrams are dis- 
cussed in Sect. 4 and our conclusions can be found in Sect. 5. 



2. General considerations on the analysis of the 
O-C diagrams 

Finding a common explanation of the possible complex be- 
haviour of O-C diagrams is not an easy task as it requires 
careful detailed analysis and modelling of the data set. It is 
well-known that the quite general model of O-C data as the 
modulation of a secular (parabohc) trend with some periodic 
additional variations due to different physical effects can usu- 
ally describe most of the observed O-C diagrams. In this repre- 
sentation, the secular part is connected to long-term evolution 
effects (mainly to mass exchange between the components) and 
the periodic modulations can be connected to e. g. LITE caused 
by a distant third companion or, if no reasonable third-orbit so- 
lution can be given, then some quasi-periodic magnetic effects 
or repeated mass-exchange processes. 

Nevertheless, there are only a few systems where the pres- 
ence of a third companion has already been proved by inde- 
pendent measurements. In such systems (e. g. VW Cephei, 
XY Leonis), the tertiary components have a period several 
times shorter than the total observational interval of the eclips- 
ing minima range, and these are clear cases of straightforward 
interpretation. Difficulty arises if the period of the cyclical term 
(if any) is close to or longer than the observational window, 
especially if any secular behaviour is also present. It is evi- 
dent that an independent fitting of the secular parabohc term 
and sinusoidal terms - as was done in several earUer works - 
may give false results, and even a simultaneous fitting does not 
necessarily improve the situation. Instead of the detailed de- 
scription of these problems, which are discussed and pubUshed 
elsewhere, we here give only two illustrations to these state- 
ments. We show that such a simultaneous fit in the case what 
is of already one simple abrupt period change can produce a 
reasonable - secular plus periodic - orbital period change so- 
lution. 

In order to demonstrate this, we constructed different arti- 
ficial times of minima data series. All of them have the follow- 
ing ephemeris: MINi = 2445 000.0 + 0.4556788£ for E < 0, 
and MINi = 2445 000.0 + 0.4556790£' for E > 0. The data 
were generated with a random accuracy of about 0.0008 days 
and with random epochs (i. e. random E value). More accurate 
points than 0.0002 days were weighted as photoelectric data. 
Here we present two data series, one which was calculated on 
the synmietiic interval -80 000 < £ < 80 000, while the sec- 
ond one on the asymmetric one of -1 10 000 < £ < 50 000. 
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Fig. 1. The artificial O-C curves with simultaneously fitted least- 
square parabolas and Fourier-curves (above) and the period curves 
(bottom). The plotted time-intervals are twice longer than the lengths 
of the data series. 



The data series were analysed in a similar manner to the one 
described in Borkovits et al. (2002), i. e. after calculating the 
O-C curve with a preliminary Unear ephemeris, the final rep- 
resentation of the O-C was searched for by a weighted Unear 
least-squares fit in the form: 

f = Co + CiE + C2E^ + ^ fly sin jViE + bij cos jViE, (1) 

where < i < 5, 1 < 7 < 4. The frequencies v, were kept fixed 
during the individual LSQ runs, but an interval of frequencies 
was scanned and the best parameter set was selected according 
to the smallest test. 

In the present cases ; - 1, j = 1..3 were applied, e.g. 
one fundamental frequency; and its first two harmonics were 
used to obtain the very good fits, which are presented in Fig. 1. 
Furthermore, what is perhaps more suprising is that the ratio 
of the Fourier-coefficients, which belong to different harmon- 
ics, almost exactly satisfy the condition of a Ught-time orbit. 
Strictly speaking, we could find reasonable third body solu- 
tions in both cases, which were almost independent of wheffier 
they were calculated only from ffie a\^2, ^1,2 coefficients (Kopal 
1978, Chap. V), or the coefficients of the second harmonics 
(Borkovits & Hegediis 1996) were also included. 

Of course, the false results can be understood considering 
ffiat removing the parabolic trend from the data set the residual 
variations form a sawtooth- hke periodic function, which can be 
represented quite well by some low order terms of its Fourier- 
series. On ffie other hand, the solutions were significantly dif- 
ferent for the two curves, which sheds light on another problem 
of ffie simultaneous fitting. Namely, if only approximately one 
period is covered by the observations, ffie coefficient obtained 
for the quadratic term and ffie period of ffie periodic variation 
are strongly coupled, and also depend on the length of ffie data. 

In what follows we analyse the O-C curves of several con- 
tact binaries, some of which clearly demonstrate this latter 
statement. However, we have to emphasise ffiat in our cases 
ffiese uncertainties are connected only to the longer period vari- 
ations but not to ffie shorter ones which are the main targets of 
our discussion here. 
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Table 1. Main parameters of the investigated stars. The spectral types Table 2. LITE solution for AB And 
are taken from SIMBAD Database^. 



z^INi = 2436109.265216(2) + 0.33 1 8897(25)£ + 3.99( 1) x Ur 



Name Sp P 1 Mi M2 Ri R2 T, T2 A refK 



E (day) 22 742 



Mo Mo Ro Ro K K Rq e' 0.10(1) 

-TD' C) 315(7) 



ABAnd G5r 0.33 W 0.60 1.04 0.78 1.03 5450 5798 2.37 1 ^ 2 436 358('448) 

OOAql G5V 0.51 A 1.04 0.88 1.39 1.29 5700 5560 3.33 2 ^ ■ ' aU 3 93f41 

DKCyg A6V 0.47 A 1.74 0.53 1.71 0.99 7351 7200 3.34 1 ^r^'"' 0158(5) 

V566 0ph FAV 0.41 A 1.56 0.41 1.51 0.86 6700 6618 2.91 3 i. ° -1 

UPeg G2V 0.37 W 1.15 0.38 1.22 0.74 5860 5785 2.52 4. ^ 



1.90 



1: Baran et al. (2004); 2: Hrivnak (1989); 3: Niarchos et al. (1993); 4: 
Pribulla & Vaiiko (2002) 

a: More recently Pych et al. (2004) found a mean spectral type of G8V 



-W3(Mo) i' = 90° 0.40 
i' = 60'' 0.48 
i' = 30'' 0.95 



3. Analysis of real systems 

We concentrated mainly on systems whose observation covers 
more than a half, or even a complete century. The most im- 
portant parameters of the analysed five systems are listed in 
Table 1 . The minima list of the five stars are not included in 
this paper due to their enormous size. Nevertheless, the total 
collection can be requested from the authors. 

In several cases earlier observations are photographic or, 
more frequently, less accurate patrol measurements or visual 
observations. They were taken into account as well, at least at 
the first step of our analysis, but with different weights. We used 
four different weights: 1: visually observed minima; 2: plate 
minima; 5: photographic normal minima; and 10: photoelectric 
observations (both photomultiplier tube and CCD). 

3.1. AB Andromedae 

This overcontact binary that belongs to the W subclass has been 
very well observed since its discovery almost 80 years ago. A 
comprehensive analysis of the photometric data obtained be- 
tween 1968 and 1995 was recently published by Djurasevic et 
al. (2000), who concluded that the light-curve variations can be 
explained well by variation of the spot activity. Two detailed 
period investigations were carried out for this system in the last 
decade. Although two different methods were applied in those 
papers. Both Kalimeris et al. (1994) and Borkovits & Hegediis 
(1996) detected the combination of secular period variation and 
cyclic period change with a period which is close to the length 
of the observational interval. This fact makes it especially in- 
teresting to carry out a new 0-C study with an extended data 
series which is longer by about 10 years. 

The analysed O-C data now cover 37 181 days from HJD 
2416 103 to HJD 2453 284. As the photographic and photo- 
electric data cover all the interval well, apart from the approx. 
6 000 day-long gap after the first three data, we omitted the vi- 
sual data in our analysis. A weighted LSQ fit of Eq. 1, with one 
fundamental frequency and with its first harmonic, resulted in a 
very reasonable fit. Nevertheless, after subtraction of this solu- 
tion we found a second, small-amplitude quasi-periodic mod- 
ulation with a period of Pmod ~ 7 000 days. In order to inves- 
tigate this possible modulation, we repeated our LSQ fitting 



0.60 




-0.40 ^ 

10000 20000 30000 40000 50000 60000 70000 

JD - 2400000 

Fig. 2. The O-C curve of the eclipsing binary AB And, with the lin- 
ear ephemeris taken from the GCVS (Kholopov et al. 1998) (above), 
and the simultaneous fit with two fundamental frequencies. We plotted 
the visual points, too, although these were discarded from the fitting 
procedure (see text for details). 

procedure with two fundamental frequencies. At the smaller 
frequency (i.e. longer period) we also took the first harmonic 
into account, while at the higher frequency we did not fit any 
harmonics. Our final solution can be seen in Fig. 2. Supposing 
that the longer period periodic term arises from a light-time or- 
bit, the orbital elements of the wide orbit of the binary around 
the centre of mass of the hierarchical triple system can be cal- 
culated easily. The results are listed in Table 2. Nevertheless, 
the calculated minimal mass of the hypothetical tertiary sug- 
gests that, if the tertiary is not a degenerate object, it has to 
produce a significant amount of third light, which should be 
observed. In contradiction to this statement, none of the pre- 
vious light-curve solutions (at least those known for us) con- 
tains any third Ught. On the other hand, we have to keep in 
mind that, due to the complex effect of the individual fitting 
parameters on the light curve solutions of contact eclipsing bi- 
naries, one can obtain similarly satisfactory approximation of 
the same problem both with and without assuming third light; 
thus it would be necessary to reanalyse the Ught curves of AB 
And with third Ught. Furthermore, slight indirect evidence for 
the motion of the binary in a triple system arises from the com- 
parison of the systemic velocity (Vo) of the system in the lat- 
est two radial velocity surveys. While Hrivnak (1988) obtained 
Vo = -24.6 ± 0.9 km s"' from the JD 2 445 886-2 446285 mea- 
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Table 3. LITE solution for OO Aql 



1940 1960 
Years 




15000 20000 25000 30000 35000 40000 45000 50000 55000 
JD - 2400000 

Fig. 3. Top: The LITE fit with the period P' = 22742 days. In order 
to see it better we removed the parabolic term, as well as the shorter 
period sinusoidal component. Bottom: The short term (/"mod = 6 695 
days) modulation on the O-C curve. For the better seeing we removed 
the other components of our final solution. 



2.50 




15000 20000 25000 30000 35000 40000 45000 50000 55000 60000 
JD - 2400000 

Fig. 4. The calculated systemic radial velocity (Vo) variation of AB 

And according to the LITE solution listed in Table 2. We marked the 
positions of the radial velocity measurements of Hrivnak (1988) and 
Pych et al. (2004). See text for details. 



surements, the most recent study of Pych et al. (2004) around 
JD 2 452 500 resulted in Vq = -27.5 ± 0.7 km s"' . In Fig. 4 we 
plotted the calculated systemic radial velocity variation of the 
system due to the revolution around the centre of mass of the 
triple system. As can be seen, our solution suggests a AVq » -1 
km s"^ variation between the above-mentioned two sets of ob- 
servations, which is not all that far from the measured differ- 
ence. Another possibility is to explain the cycUc variation by 
magnetic activity (see Sect. 4.3). 

3.2. 00 Aquilae 

This high mass-ratio A-type overcontact binary was discov- 
ered by Hoffleit (1932). It belongs to the small group of longer 
period contact binaries which are supposed to have evolved 
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24438613.18013(1) + 0.5067920(14)£; 
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into the contact phase in the (astronomically) near past (see 
e. g. Hrivnak et al. 2001, and references therein). Spectroscopic 
study by Hrivnak et al. (2001) revealed strong Mg II h & k 
emission lines varying with time, as clear evidence of chromo- 
spheric activity, although they also pointed out that OO Aql is 
a less active system than other W UMa-type ones. The most 
detailed period study in the recent literature was published by 
Demircan & Giirol (1996). They concluded that the most prob- 
able two scenarios for the large structure period variation were 
an abrupt period jump with AP » -0.6s around 1963 or a sinu- 
soidal variation with a period of P' ~ 90 yr, although this sec- 
ond statement was quite ambiguous, as the total curve covered 
only one maximum without any minima. They also mentioned 
the possibility of a parabolic representation. Nevertheless, in 
the present evolutionary state of the system, it could be ex- 
pected that the binary evolves toward a smaller mass-ratio state, 
which should have resulted in an increasing period as was 
clearly in contradiction with the 1995 state. Demircan & Giirol 
(1996) also investigated the presence of shorter, small ampli- 
tude fluctuations and concluded that these variations did not 
show simple periodic nature. 

Since this aforementioned study, the picture has changed 
significantly. Around 1996 the period of the system also varied 
(see top panel in Fig. 5). Our O-C data extend on a 25 919 day- 
long interval between HJD 2 426 892 and 2 452 8 1 1 . Although 
it cannot be excluded that this is a second abrupt period change, 
we may also suspect that the large amplitude sinusoidal vari- 
ation reached its minimum. Consequently, we applied our fit 
with two fundamental periods: one (with its first harmonic) for 
the long period sinusoidal variation, and one for the short pe- 
riod fluctuations, which clearly reveals periodic behaviour in 
contrast to the statement of the previous authors. The simulta- 
neous fit is plotted in the top panel of Fig. 5, while in the bottom 
panel we iUustrate the periodic nature of O-C lacking long pe- 
riod terms. Assuming that the P' « 75yr variation arises from 
the presence of a more distant tertiary component, the orbital 
elements of the binary in the triple system are Usted in Table 3. 
As seen there, the minimal mass of the tertiary is again so large 
that (if not a degenerate object) it should be seen via the light 
curve with significant third Ught contribution, as well as via 
its spectral lines. Nevertheless, the previous Ught-curve studies 
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Fig. 5. The O-C curve of the eclipsing binary OO Aql with the calcu- 
lated linear ephemeris and the simultaneous fit with two fundamental 
frequencies (above), and the smaller period fluctuations after the sub- 
traction of the longer period curve. The Fourier-fits were carried out 
simultaneously. 

give clear indirect evidence against the presence of any sig- 
nificant third light; namely, that the inclination of the system 
without the supposition of any third light was found to be close 
to i = 90° (Hrivnak 1989; Djurasevic & Erkapic 1999). As the 
third light reduces the depth of the light minima, a given depth 
could be reproduced in its presence only by a larger inclination 
value than in the case of no third light, which would be physi- 
cally unrealistic for the present star. Consequently, even if the 
75 year-long periodic amplitude is real, an alternate explana- 
tion is needed for it. 



3.3. DK Cygni 

Although the discovery of this A-type overcontact binary was 
reported in the same paper (Guthnick & Prager 1927) than as 
was the above-mentioned AB And, and their brightness is also 
similar, this system was less popular with the observers. Paparo 
et al. (1985) found a continuous period increase in the system, 
which was confirmed by Wolf et al. (2000). No spot activ- 
ity was observed in this system as a survey of the literature 
showed. 

The collected times of minima data in our analysis cover 
28 542 days from HJD 2 424 760 to HJD 2 453 302. In the anal- 
ysis we also used seven recently obtained unpublished minima, 
listed in Table 4. Omitting the visual data, a weighted LSQ fit 
resulted in the following parabolic ephemeris (see Fig. 6): 

MINi = 2451000.1031-i-0.470693909£ + 5.862xl0-'i£2 (2) 

which is in good agreement with the previous results. Despite 
this fact, we think that it is not without benefit to reanalyse 
the O-C of this system, as well as any other eclipsing bina- 
ries, every 5-10 years by the use of the most recent data points. 
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Fig. 6. The O-C curve of the eclipsing binary DK Cyg, with the linear 
ephemeris taken from Kiss et al. (1999) (above), and the residual curve 
after subtracting the parabola. 

Table 4. New times of minima of DK Cyg observed by R. Csorvasi 
with CCD cameras mounted on the 40-cm Cassegrain telescope at 
Szeged Observatory, and the 90-cm Schmidt telescope at Piszkesteto 
Mountain Station of Konkoly Observatory. 



HJD 



Error Type 



2453223.4286 ±0.0011 s 
2453228.3681 ±0.0008 p 
2453 246.4950 ±0.0013 s 
2453 247.4346 ±0.0007 s 
2453285.326 ±0.002 p 
2453 286.2657 ±0.0007 p 
2453 302.2672 ±0.001 p 



This is first because, in several cases, it cannot be decided from 
observations with gaps whether an abrupt period jump has oc- 
curred in a given time or if we see some continuous variation. 
However, a more or less continuous monitoring could soon dis- 
cover if something interesting is happening, and would lead 
the observers (nowadays sometimes skilled amateurs with ad- 
vanced equipment like CCD cameras) to observe the given bi- 
nary more frequently. Second, returning to this particular sys- 
tem, the constancy of the period variation (at least in the last 80 
years), and the lack of any other features from the O-C curve 
needs some explanation, especially because in most cases the 
O-C curves of the contact binaries show more complex be- 
haviour. 

3.4. V566 Ophiuchi 

This A-type contact binary was discovered by HofFmeister 
(1935). Although it was one of the most popular contact bi- 
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Table 5. LITE solution for V566 Oph 



MINi = 2440418.49522(7) + 0.4096433(22)£; + 1.55(11) x IQ-^^E^ 



P' (day) 
e' 

OJ' (°) 

t' {HJD) 

an sin i' AU 

/(ms) Mo 

Ki2 kms"' 



7 250 

0.31(8) 

180(14) 

2453 590(267) 

0.63(5) 

0.00065(13) 

1.00 



m3(Mo) 



= 90° 0.14 
= 60° 0.17 
= 30° 0.30 



naries during its "first 60 years", it seems to have been ne- 
glected in the past decade. V566 Oph did have a stable light 
curve and constant period in the past (Bookmyer 1976; Eaton 
2001). After 1983 the depths of eclipse were increased (Lafta & 
Grainger 1985). The behaviour of V566 Oph is not understood, 
so it requires further study. Due to the light curve changes, the 
system seems to be a rather active star 

The latest period analysis was carried out by Qian (2001b), 
who fitted a least- squares paraboUc ephemeris and, after sub- 
tracting this parabola, he also fitted a pure sinusoidal term, with 
an amplitude of O?0037 d and a period of 20.4 yr. He interpreted 
this latter component as the trace of a possible third companion. 

As already mentioned in Sect. 2, the non-simultaneous fit- 
ting of the secular and periodic terms might give inaccurate 
results. This was the main reason we repeated this earlier O-C 
analysis; the other reason, of course, is that some new minima 
times have been obtained since then. Our times of minima data 
cover the interval HJD 2 434 179-2 452 854. As was mentioned 
earlier, there are only very few recent minima. We only found 
two within the last 7 years. In the present analysis we omit- 
ted the visual minima, because first the total observing interval 
is satisfactorily covered by photoelectric minima and, second, 
the number of the visual minima is so small that their statistical 
significance, which seems to be evident e. g. in the case of AB 
And or U Peg, fails here. 

After omission of four evidently bad photoelectric data, we 
simultaneously fitted a quadratic polynomial and the Fourier- 
terms up to the first harmonics. The best fit is plotted in Fig. 7, 
while the calculated LITE orbit parameters are hsted in Table 5. 
Comparing our results with those of Qian (2001b), both the 
parameters of the parabola and the periodic terms are simi- 
lar, for two reasons. First, the period of the quasi-sinusoidal 
term is less than half of the length of the data series; conse- 
quently, the periodicity can be calculated almost independently 
of the quadratic term parameters. Second, the amplitude of this 
term is small with respect to the secular period change rate. 
Although in Table 5 we present the LITE solution for this al- 
most 20-year period modulation, a possible magnetic origin is 
also investigated in Sect. 4.3. 
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Fig. 7. The O-C curve of the eclipsing binary V566 Oph, with the lin- 
ear ephemeris taken from the GCVS (Kholopov et al. 1998) (above), 
and the LITE fit on the parabola-subtracted curve. The quadratic and 
Fourier-fits were carried out simultaneously. 



3.5. U Pegasi 

This W subtype system is one of the oldest known W UMa type 
binaries, as was already discovered 110 years ago (Chandler 
1895). In a recent study, Djurasevic et al. (2001) investigated 
the long term photometric variations of the system over the 
time interval of 1950 to 1989. They concluded that the complex 
variations of the light-curves during these four decades can be 
explained by the variable spot activity of the cooler component. 
The latest period study was carried out by Pribulla & Vanko 
(2002), who analysed the features of the O-C diagram between 
1894 and 2001. These authors explained the century-long pe- 
riod decrease by conservative mass transfer from the more mas- 
sive to the less massive component. Furthermore, they found a 
short period (Pmod - 18.8 + 3.9 years), low amplitude oscilla- 
tion. This feature had also been reported earher by Zhai et al. 
(1984). 

Our O-C data cover 39 856 days between HJD 2413 094 
(1894) and 2452950 (2003). These data are plotted in the up- 
per panel of Fig. 8. We carried out our analysis in a similar 
way as the one described in the case of AB And; i. e. we fitted 
simultaneously a quadratic polynomial and two fundamental 
frequencies allowing the first harmonic at the longer period as 
well. Our best fit gave the following two periods: Pi » 31 059 
days, P2 ~ 6 544 days. If the first period belongs to a light-time 
orbit due to the existence of an additional star in a hierarchical 
triple (multiple) system the corresponding orbital elements and 
derived quantities can be read in Table 6. The latter period is 
approximately the same as that reported in the previous papers; 
what perhaps is more interesting, this is also very similar to the 
length of the modulations in the previous stars. Nevertheless, in 
the present case it is questionable that this is a real period vari- 
ation. As can be seen well in Fig. 9, the total amplitude of this 
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Table 6. LITE solution for U Pes 



MINi = 2436511.6698(5) + 0.3747809(5)£: - 3.61(10) x 10-'i£? 



P' 
e' 

OJ' 

t' 

Ugu sin /' 
/(ma) 

K\2 



(day) 

n 

(HJD) 
AU 
Mo 
kms"' 



31059 
0.38(5) 
112(9) 

2445 935(569) 
1.96(13) 
0.00105(20) 
0.74 



m3(Mo) = 90° 0.14 
i' = 60° 0.17 
!' = 30° 0.30 
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Fig. 9. The shorter period sinusoidal fit on the 0-C of U Peg. For the 
sake of clarity, we subtracted the parabola, as well as the longer period 
LITE solution, from the original curve. 




Table 7. Derived parameters from the period variations 



10000 15000 20000 25000 30000 35000 40000 45000 50000 55000 

JD - 2400000 

Fig. 8. The O-C curve of the eclipsing binary U Peg with our final 
solution, plotted against the linear ephemeris taken from the GCVS 
(Kholopov et al. 1998) (above), and the pure LITE fit on the parabola 
subtracted curve. The quadratic and Fourier-fits were carried out si- 
multaneously. 

"variation" is only A2 ~ 0^003, i. e. smaller than the scatter of 
the photoelectric minima, which is usual in W UMa binaries, 
as a consequence of the varying light-curve shapes distorted by 
spots (see e. g. Kalimeris et al. 2002). In our opinion it might 
be imagined that the presence of such a periodicity, especially 
between 1960 and 1980, is only a result of the small number of 
photoelectric measurements. However, the fact that similar pe- 
riodicity was found in the other cases, too, makes it reasonable 
that this feature might be real. 

4. Discussion 

The five investigated systems are distributed in the range of ef- 
fective temperature from T^g = 7400 down to T^s = 5400, 
which is the transition region where a convection zone devel- 
ops in the envelope of the stars toward the lower temperatures. 



Name 




\m\ 


^ mod 






xlO-' (y- 


') xlO-^(Moy- 


') (d) 


(d) 


AB And 


2.65 


1.25 


6695 


0.0020 


00 Aql 






7467 


0.0036 


DK Cyg 


1.93 


0.49 






V566 0ph 6.76 


1.25 


7 250 


0.0035 


UPeg 


-1.88 


0.36 


6544 


0.0016 



Analysis of the O-C diagrams of the systems has revealed dif- 
ferent periodic variations in four cases besides the secular de- 
velopment of the orbital period. 

4.1. Secular period changes 

Four of the five O-C curves show evidence of continuous or- 
bital period change with a constant rate during the total ob- 
serving interval. The only exception was OO Aql, although it 
is possible that in this system only one (or two) abrupt period 
jump(s) obscures this behaviour. 

The most usual explanation of this kind of secular period 
variation is mass exchange in the system. The widely used ap- 
proximating formula for the calculation of the mass exchange 
rate is as follows: 



-OT12 



<1 



l-q^3P 



-ntn 



q 2c2 



(3) 



where ci, C2 directly come from the O-C ephemerides in the 
form of Eq. (1). Note that the equation requires constancy of 
the total mass, as well as of the total angular momentum. The 
rate of the secular period change, as well as the calculated mass 
exchange rate, are tabulated in Table 7. 

Nevertheless, it is necessary to note, that in the case of 
AB And and U Peg, the coefficients of the quadratic terms and 
the period and amplitude of the longer period quasi-sinusoidal 
period variations are not independent of each other as we will 
see later. 
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4.2. Longer period cyclic variations 

In three cases the simultaneous fitting gave longer, as well as 
shorter, period cycles. These longer scale variations are roughly 
62, 75, and 85 years long and naturally could be identified as 
a LITE, i. e. effect due to the presence of a third body. The 
corresponding LITE solution parameters were given in the cor- 
responding tables along with the detailed discussion of the sys- 
tems. The identification of these periods as LITE solution can 
be considered only as a simple designation, and their real origin 
should still require a careful analysis of all cases. 

It should also be emphasised that these periods are in the 
same order as the total interval of the observations, and the sep- 
aration of the different effects in the analysis are doubtful. For 
example, let us consider e. g. the case of AB And. In this case 
the fundamental Fourier-terms are as follows: 

/ = 0.015 sin(9.17 x lO^^E) - 0.017 cos(9.17 x IQ-^E) + ...,(4) 

which can be written as 



/ = -0.017 + 1.38 X IQ-^E + 7.15 x lO""^^ -i- 



(5) 



where as seen well, the quadratic term is in the same order of 
magnitude than the coefficient of the parabolic one. This fact 
is critical mainly in these two abovementioned cases (AB Cas 
and U Peg), as in both systems only the last half quasi-sines are 
covered well by good-quality photoelectric data. On the other 
hand, despite this fact, we can stress that at least the order of 
magnitude of the calculated secular period change is believed 
to be correct. 

4.3. Shorter period cyclic variations - manifestation of 
magnetic activity cycles? 

It is an interesting fact that four of our five systems show small 
ampUtude cycUc fluctuations with very similar periods (see 
Cols. 4-5 in Table 7). The only exception is DK Cyg, the hottest 
in the sample, with A8 spectral-type components, whilst the 
other four are later than FO down to G5. 

These modulation periods range from 6 500 to 7 500 days 
(18-20 years), which are quite similar to those observed in 
other types of stars showing magnetic activity. Comparing 
them to the longer period components of the detected varia- 
tions, where the observations cover only partly a cycle, these 
shorter cycles are well-covered (several times, 2-4 cycles), as 
can be seen from the discussion of the individual systems. 
Hence, the existence and characteristics of these cyclical vari- 
ations in the orbital period are well proved independently from 
any doubts about the other - longer period - components. 

The startUng resemblance of orbital period variations in 
these more or less similar systems suggests some parallel in 
their origins and makes some common physical explanations 
Ukely. In our opinion this feature might be indirect evidence 
of similar magnetic cycles in the investigated binaries, and the 
detected variations might be caused by it. 

We mainly refer to Lanza & Rodono (1999), who found a 
relation between the orbital period of a binary system and the 



magnetic activity cycle by supposing synchronisation between 
the orbital motion and rotation of the member stars: 

log^mod [yr] = 0.018 - 0.36(±0.10)log|^[sec]. (6) 

^orb 

For the binaries studied here, this formula predicts f mod ~ 
7 900, 9 200, 8500, and 8 200 days, respectively (note values 
very close to the length of the Sun's magnetic cycle). These 
values are on the same order of magnitude as the observed 
ones. This supports our conjecture that magnetic activity cy- 
cles could explain the observed short term orbital period varia- 
tions. Furthermore, beyond our sample, Awadalla et al. (2004) 
reported a similar periodicity of Pmod ~ 6 500 days in the O- 
C curve of the late-type (F8) W UMa binary AK Her, giving 
additional support to the explanation above. 

The rate of this period variation can be expressed easily 
(see e. g. Applegate 1992) as 



AP - Ao-c 

" ^mod 



(7) 



where AP/P is the relative period variation, Aq c the half- 
amplitude of the O-C variation, and /"mod the cycle length of 
the magnetic activity. 

In the case of AB And, applying f mod = 6 700 days and 
A = 0.002 days, respectively, we get AP/P = 1.9 x 10 ''d/d. 
The variation of the gravitational quadrupole momentum is 
(Applegate 1992): 



1 , . AP 

AQ = --MRHR/a)-^ — 



(8) 



where M is the mass of the active star, R the radius of the active 
star, and a = 2A the separation between the components. A is 
the semi-major axis, and note that in contact binaries the or- 
bit is circular. A calculation with the parameters of the primary 
and the secondary (see Table 1) gives AQi = 4.02 X lO^^gm^ 
and AQ2 - 6.97 x lO^^kgm^, respectively. These values cor- 
respond to the typical values in active binary stars (Lanza & 
Rodono 1999). Due to the similar absolute dimensions of the 
other investigated binaries, the same calculations give physi- 
cally realistic results for them as well. 

Nevertheless, we also have to emphasise that in a recent 
study Qian (2003) did not find clear evidence of decades time- 
scale period variations in similar systems. He investigated the 
period variation of five contact binaries (GW Cep (G3), VY Cet 
(G5V), V700 Cyg (G2V), EM Lac (G8V), AW Vir (GOV)) and 
found a secular period change in every case, and he discovered 
that only two of them show cyclic period variations (VY Cet 
(P3 = 7.3 years), V700 Cyg (P3 = 39.8 years)). It is an inter- 
esting fact that only 40 percent of these G spectral type sys- 
tems have cyclic period variations, despite expecting such a 
phenomenon in all cases. Although both samples (that of Qian 
2003 and our) are small, and therefore not representative, the 
controversion should be studied. Our own opinion is as follows. 
The O-C diagrams of the three remaining systems (GW Cep, 
EM Lac and AW Vir) consist in only a few points, and there are 
large gaps between the data. This is due to the fact that these 
systems were underobserved in the past; therefore, the present 
observational material is enough to detect the secular period 
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variation only - in spite of applying all available minima. Gaps 
and scattered observations (Uke visual and photographic ones) 
obstruct the demonstration of such cyclic variations. 

As mentioned already, the spectral type of the systems is 
later than FO, where we have found short period variations in 
the 0-C diagrams, the earlier type DK Cyg does not show this 
feature. This type of behaviour is pretty similar to the observed 
one among the Algol-type RS CVn systems (see Hall 1989) so 
that the magnetic activity can be observed only in systems later 
than F5. Because the magnetic dynamo theories are based on 
a strong connection with the presence of a convection zone in 
the outer envelope predicting magnetic activity for them. Since 
the convection zones in the envelope appear in lower mass MS 
and cooler stars, this type of relation is a trivial requirement 
but, of course, the strict border lines can be different. Although 
our sample is not a systematic one, it indicates the existence of 
a separation border somewhere around the spectral type FO. 

5. Conclusions 

The orbital period is very sensitive to any effect and it could 
be measured with very high precision down to 10"^ part of the 
period, and we have a lot of series of minimum observations, in 
some cases back to the 19th century. In this study we have in- 
vestigated the 0-C diagrams of five W UMa systems from A8 
to G5 spectral type components. In four cases a short time scale 

20 years) cyclic variation of the orbital period was found. 
All these systems contain later than FO spectral-type compo- 
nents, which means that they have deep convection zones. As 
we showed the very similar modulation periods that were ob- 
served seem to be independent of any LITE caused by a third 
body, but a consistent explanation of their cause might also be 
a magnetic activity cycle. 

The mechanism proposed by Applegate (1992) seems to be 
responsible for the observed variations, and the investigation 
of close binary systems from this aspect could be an important 
tool for studying the magnetic properties of stars in general. 
However, contact binary stars are generally faint and the ge- 
ometry very complex, so that it is very difficult to apply these 
methods until further further developments (see e. g. Hendry & 
Mochnacki 2000). 

Of course, these studies have to be complemented by de- 
tailed analyses of the systems from other viewpoints, e. g. ob- 
servations of the time variable features of their spectra, mod- 
elhng of time development of their light curves, etc. For ex- 
ample, to check our hypothesis about the magnetic activity, we 
should follow in time the spotted hght curve models of the sys- 
tems determining their spottedness-index and its correlation 
with the observed 0-C variations. These works are in progress 
and will be pubhshed in due time. 
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